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Summary. — We study the radiation damping effect on the relativistic acceleration 
of electron-positron plasmas with two-and-half-dimensional particle-in-cell (PIC) 
simulation. Particles are accelerated by Poynting flux via the diamagnetic rela- 
tivistic pulse accelerator (DRPA), and decelerated by the self-consistently solved 
radiation damping force. With £l ce /uj pe > 10, the Lorentz factor of the highest 
energy particles reaches 7 > 100, and the acceleration still continues. The emit- 
ted radiation is peaked within few degrees from the direction of Poynting flux and 
strongly linearly polarized, which may be detectable in 7-ray burst(GRB) observa- 
tions. We also show that the DRPA is insensitive to the initial supporting currents. 
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One of the unsolved problems in astrophysics is the acceleration of nonthermal high- 
energy particles, whose radiation is observed from pulsars, blazars, gamma-ray bursts 
and black holes. Recently, a new mechanism of relativistic nonthermal particle accelera- 
tion, called the Diamagnetic Relativistic Pulse Accelerator (DRPA), was discovered using 
multi-dimensional Particle-in-Cell(PIC) simulations. When a plasma-loaded electromag- 
netic pulse expands relativistically, the self-induced drift current creates ponderomotive 
trap, which drags only the fast particles in the trap and leave slow ones behind. At 
late times the DRPA reproduces many of the unique signatures of GRBs. including time 
profiles, spectra and spectral evolution^. 

When charged particles suffer extreme acceleration, radiation loss and damping can 
become important in the plasma energetics and dynamics. However, conventional Particle- 
in-Cell (PIC) simulations of collisionless plasmas have not included radiation effects. In 
this article we report PIC simulation results using a newly developed 2-1/2-D code that 
includes self-consistent radiation damping. 

In PIC simulations, it is impractical to include high-frequency radiation into the 
electromagnetic field calculation because the radiation wavelength is much shorter than 
the spatial resolution of the fields (<~ Debye length Arj = c/oj pe , where ui pe = y/ '4ir pe / 'm e 
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is the electron plasma frequency). Accelerated particles can emit up to the critical 
frequency lu c — 3^/ 2 Q cei where 7 = E/m e c 2 = 1/y/T— v 2 /c 2 and VL ce = eB/(m e c) is 
the electron gyro-frequency. The ratio of the critical radiation wavelength A c to A^ is 
given by A c /Ad = (27nj pe )/(37 3 r2 ce ), which is <gC 1 because uj pe /Q, ce < 0.1 in magnetic- 
dominated cases and 7> 1. 

Instead of calculating the high-frequency component directly, we introduce a radiation 
damping force in the form of the Dirac-Lorentz equation The relativistic damping 
force term f rad is given by (see 3^ for detailed calculations) 

f ,L 

J rad ~ ^rad S I 

(1) + -Bx(Bxv) + -E(v ■ E) 

c c 

where v is the velocity, and E and B are the self-consistent electric and magnetic fields. 
Here we introduce a non-dimensional factor k ra d given by 

(2) k rad = = 1.64 x 1(T 16 x S(gauss), 

c 

where r e — e 2 / (mc 2 ) is the classical electron radius. The first square bracket term of the 
radiation damping force (JJJ represents the radiation damping due to the ponderomotive 
force acceleration. The third square bracket term is Compton scattering by large scale 
(A > Ad) electromagnetic field which reduces to Thomson scattering in the classical 
limit [5]- We should note here that the scattering between high frequency radiation and 
particles is not considered since all fields are averaged over Debye length. 

The power spectrum analysis shows that our simulation model with ui pe /fl ce = 10 
corresponds to the GRBs with initial magnetic energy ~ 10 51 ergs by assuming a 4ir shell 
of thickness 10 12 cm and radius 10 13 cm. If GRBs are originated from a region < 10 8 cm, 
this magnetic energy implies an initial B > 10 13 G, which corresponds to k rad ~ 10~ 3 . 

We should restrict ourselves not to reach the quantum-limit, Ml ce ~ m e c 2 or B > 
4.4 x 10 13 G, which corresponds to k ra( i = 7.2 x 10~ 3 , or the formula QJ fails. We 
choose k ra d from zero to 10~ 3 in the simulation to enhance the radiation effect and 
\frad\ T sim — | F ext | ^ ce s0 we can see * ne difference between radiative and non-radiative 
(krad = 0) case within the simulation time-scale T S i m = O(10 4 ^l~^). 

We use the 2-1/2D explicit PIC simulation scheme using the explicit leap-frogging 
method for time advancing [Hj. Spatial grids for the fields are uniform in both x and z 
directions, Ax = Az = \r>. The simulation domain in the x — z plane is —L x /2 < x < 
L x /2 and < z < L z with a doubly periodic boundary condition in both directions. 

Following Liang et. al. PQ, the initial plasma is uniformly distributed at the center 
of the simulation box, — 6Ax < x < 6 Ax and < z < L z . The background uni- 
form magnetic field Bq = (0,i?o,0) is applied only in the same region, so that the 
magnetic field freely expands toward the vacuum regions, x > 6Ax and x < — 6Ax 
with accelerating plasma. We choose L x to be long enough so that plasma and EM 
wave never hit the boundaries in the x direction within the simulation time. The ini- 
tial temperature of plasma is assumed to be a spatially uniform relativistic Maxwellian, 
kBT e = ksTp = lMeV, where the subscripts e and p refer to electrons and positrons. 
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Fig. 1. - An example ol the DRPA with initial current J = V x B , u) pe /Q, ce = 0.1, k ra d = 0. 
The initial magnetic field and current distribution are shown in panels (a) and (b). The phase 
plot of particles at t£l ce — 4800 is shown in panel (c). 



Since its discovery, questions have been asked if DRPA would persist if the initial 
static B field is confined by supporting currents (J = V x B), instead of being simply 
superposition of opposite traveling waves with initial ,7 = 0. Figure ^ shows an example 
of a PIC simulation with initial current distribution J = V x B. We use a finer grid 
Ax = Ad/3 in this case to resolve the gradient of fields near the plasma surface. At t = 0, 
magnetic field is smoothly decreased toward the edge x/Xo = ±6, and initial current 
is locally distributed near the edges. However, the initial current is rapidly dissipated 
by the expanding plasma, and two EM pulses are formed as in the case without initial 
current. The asymptotic phase plot is shown in Fig.H(c). It is basically identical to the 
DRPA phase plot without initial J. 

We choose six different sets of parameters shown in Table [i] for the radiative calcula- 
tions, by changing k ra d = 0, 10~ 4 , 10~ 3 and oj pe /£l ce = 0.1, 0.01, and run simulations 
for each case. 

Figurc[5]shows the momentum distribution of particles for (a) Run A(x < 0) and Run 
C(x > 0), and (b) Run D(x < 0) and Run F(x > 0). The DRPA accelerates electrons 
and positrons in the same direction along the x axis, whereas electric field accelerates 
electrons and positrons oppositely along the z axis, forming X shape distribution in the 
Px — Pz plane as a result. Resulted induced current J 2 accelerates particles in the x 
direction by the ponderomotive force J x B. The ponderomotive force creates successive 
'potential wells' in the x direction, which captures and accelerates co-moving particles. 
We emphasize here that there is no charge separation in the x direction because there is 



Table I. - The parameters for 


radiative runs 
k 


LOpe 1 i^ce 


Duration tQ. C e 


Run A 





0.1 


10000 


Run B 


10~ 4 


0.1 


10000 


Run C 


10 -3 


0.1 


10000 


Run D 





0.01 


70000 


Run E 


10" 4 


0.01 


70000 


Run F 


10~ 3 


0.01 


70000 
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Fig. 2. - Momentum distribution of particles for Run A [(a), x < 0] and Run C [(a), x > 0] at 
tQ. C e — 5000 (circle, asterisk) and tO ce = 10000 (cross, plus), and Run .D [(b), x < 0] and Run 
F [(b), x > 0] at tQ. ce — 35000 (circle, asterisk) and tO ce ~ 70000 (cross, plus). Results in the 
positive and negative x directions are identical in all cases. 



no mass difference between electron and positron. 

Obviously, particle momenta in both x and z directions are radiated away in both weak 
and strong magnetic field RD cases. For high energy (7 ^> 1) particles, the Compton 
scattering [the third term in Eq. becomes dominant damping force, and makes the 
DRPA acceleration less efficient. 

With radiation damping, we can calculate the self-consistent radiation field and its 
angular dependence directly from the velocity and acceleration of each particle. Intensity 
/ and polarization II of the radiation received by the observer located at x are given by 

El El 



ret 



(3) /(n,r)=£[|f 

i 

where 



, n(n,r) 



(4) 



U(n,T) = 2^2[(E i -y)(E i -z)] x 



et • 



[Elf + (ElY - 2E*E* + m 



ret , EUn,r) = ^[\E l -z 



ret 

1 

e nx[(n- f3 t ) x /3 Z ] 
c {l-n-(3fR ' 



where the sum is over all particles along each light cone, n is a unit vector in the 
direction of x — r(r), (3 = v(t)/c, and (3 = df3/dt. We assume that |a;| 3> \r\ so that 
n is parallel to x. The square brackets with a subscript "ret" mean that the quantity 
in the brackets is evaluated at the retarded time r = t — R/c, where R = \x — r\. To 
specify the direction of the observer with respect to the x axis, we introduce 9 and (j> as 
n = (cos cos 4>, cos 8 sin cj>, sin 9) . 

In Figs.|2Ia) and (b), the contour plot of instantaneous intensity before taking the 
summation over the retarded time is plotted, as a function of local time t with <f> = 9 = 0, 
illustrating the ray-tracing technique used in Eqs. J2Jl and |@}. We take a sum of intensity 
along the light cone r = t — R/c —const., which is indicated as a solid line in the panel 
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Fig. 3. - Contour plots of instantaneous intensity log 10 I per particle as a function of local time 
t for 9 = <f) — for (a) Run C and (b) Run F. Total detected intensity (solid lines, right scales) 
and polarization (dashed lines, left scales) as functions of observational time r for (c) Run C and 
(d) Run F. Intensity is in arbitrary scale. To obtain time dependence of intensity, instantaneous 
intensity from each particle is summed up along the light cone r = t — R/c =const., shown as a 
black line in panel (a). The light cone moves horizontally leftward with r. 



(a), up to t£l ce = 10000 for Run C and t£l ce = 70000 for Run F. The light cone moves up 
in t with increasing r, and we take r = when the pulse front reaches to the observer. 
We also note that the intensity is extremely asymmetric because only energetic particles 
accelerated in positive x direction can radiate strong emission to the observer. 

The time dependence of detected intensity and the polarization with <f> = 6 = are 
shown in Figs.EJc) and (d). Energetic particles are bouncing back and forth within the 
ponderomotive potential well, and slower particles are dropped off to the next well, which 
broaden the spatial distribution and the resulting radiation duration. 

Intensity / is shown as solid lines in Fig. [3Jb) and (d) , indicating the duration time 
of strong intensity radiation is 20r£! ce for Run C, and 200 ~ 400ri7 ce for Run D with 
single peak. Since the driving electromagnetic field is linearly polarized, we expect that 
the radiation is strongly polarized with the small depolarization coming from the initial 
random y velocity. Dotted lines in Fig. EJb) and (d) shows that the radiation is strongly 
linear-polarized as anticipated. 

In Fig. 01 we show the total radiation intensity J I((j>,8)dT. The radiation is a very 
short pulse in Fig. [3] because we consider single simulation box centered at the origin 
only. However, because of the periodic boundary condition in the z direction, we should 
consider the multiple simulation boxes along the z axis, and consider time delay and 
angle difference toward the observer from each simulation box, which is very complex 
even in the Cartesian coordinate. 

Instead of including these geometrical effects into the intensity calculation, we simply 
compare the total fluence integrated over i as a function of the angle. Angular fluence 
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4>, 6 (degree) 

Fig. 4. - The total radiation intensity J I(4>, 8)dr for Run C and Run F as a function of the 
angle ((f>,6). Solid: Run C, <j> = 0; Dashed: Run C, 6 = 0; Dash-dot: Run F, <f> = 0; Dotted: 
RunF, 6 = 0. 



peaks around 6* = 3 ^ 8° in cf) = cases (solid and dash-dot lines), corresponding to the 
direction of high energy particles in Fig. [2J Angular fiuence rapidly decreases with both 
<f> and 9, indicating radiation is strongly collimated in the x direction. 

Angular fiuence distribution in the 9 direction is due to the initial electric field ac- 
celeration in z-direction. Since there is no z acceleration after the EM pulse is formed, 
eventually all the z momentum will be emitted away, which narrows the intensity distri- 
bution toward positive x direction. At the same time, however, decoupling of particles 
from EM pulse in the x direction widens the distribution. Thus, the fiuence in the 9 
direction is always distributed over a finite angle. 

In summary, we observed the self-consistent radiation damping effect on the inter- 
action of the EM pulse with electron-positron plasma via a relativistic PIC simulation. 
Comparison with the non-radiative case showed that radiation damping force decelerates 
the energetic particle accelerated by the DRPA, and resulting radiation field is strongly 
linearly polarized within a finite angle from the direction of expansion both in weak and 
strong magnetic field cases, which may be detectable by 7-ray burst observations or laser 
experiments as an indication of the DRPA mechanism. 

The simulations shown here are still too short to determine the radiation pattern 
after particles are completely decoupled from EM pulse. Such questions remain as future 
problems. 

* * * 

This research is partially supported by NASA Grant No. NAG5-9223 and LLNL 
contract nos. B528326 and B541027. Authors are also grateful to ILSA, LANL, B. 
Remington and S. Wilks for useful discussions. 



REFERENCES 

[1] Liang E., and Nishimura K.,Phys. Rev. Lett, 92 (2004) 175005. 
[2] Zhidkov A. et al.,Phys. Rev. Lett, 88 (2002) 185002. 

[3] Landau L. and Lifshits E. M.,The Classical Theory of Fields, (Pegramon, New York) 
1975. 



SELF-CONSISTENT RADIATIVE EFFECT ON RELATIVISTIC ELECTROMAGNETIC PARTICLE ACCELERATION? 

[4] Jackson J. D., Classical Electrodynamics 2nd ed., (Wiley, New York) 1975. 

[5] Rybicki G. B. and Lightman A. P .,Raidiative Processes in Astrophysics, (Wiley- 
interscience, New York) 1979. 

[6] Birdsall C. K. and Langdon A. B., Plasma Physics via Computer Simulation, (McGraw- 
Hill) 1985. 



